Superconductive pumping of nanomechanical vibrations 
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Abstract 

We demonstrate that a supercurrent can pump energy from a battery that provides a voltage bias into nanomechanical 
vibrations. Using a device containing a nanowire Josephson weak link as an example we show that a nonlinear coupling 
between the supercurrent and a static external magnetic field leads to a Lorentz force that excites bending vibrations of the 
wire at resonance conditions. We also demonstrate the possibility to achieve more than one regime of stationary nonlinear 
vibrations and how to detect them via the associated dc Josephson currents and we discuss possible applications of such a 
multistable nanoelectromechanical dynamics. 

PACS numbers: 74.45.+C, 74.50.+r, 74.78. Na, 75.47.De 
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Coupling of electronic and mechanical degrees of free- 
dom on the nanometer length scale is the basic phe- 
nomenon behind the functionality of nanoelectromechan- 
ical (NEM) systems. Such a coupling can be mediated 
either by electrical charges or currents. Single-electron 
tunneling (SET) devices with movable islands or gate 



biased weak link it becomes possible to pump nanome- 
chanical vibrations in the Cooper pair tunneling region. 
As we will show below the result is a peculiar nonlinear 
NEM dynamics that affect both the supercurrent flow 
and the nanomechanical vibrations in novel ways. 



electrodes employ Coulomb forces to achieve capacitive 
0, and shuttle NEM coupling 0, Q, where the lat- 
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ter involves both capacitive forces and charge transfer. 
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Devices containing current carrying parts, on the other 
hand, will achieve NEM coupling through magnetic-field 
induced Lorentz forces. Focusing on the latter mecha- 
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nism, a simple estimate shows that for a gold nanowire 
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suspended over a few micrometer long trench, the me- 
chanical displacement due to typical currents of order 
100 nA in magnetic fields of order 0.01 T can be as large 
as one nanometer. Such displacements can crucially af- 
fect the performance of mesoscopic devices. 

In this Letter we will explore a possible scenario for 
how highly nonlinear nanoelectromechanical effects can 
arise if the magnetic-field induced electromotive force 
caused by the mechanical motion of a conducting wire 
strongly perturbs the flow of current through it. De- 
vices which contain superconductors, with their known 
extreme sensitivity to external electric fields, are the best 
candidates to achieve such strong effects and supercon- 
ducting quantum interference devices (SQUID's) that in- 
corporate a nanomechanical resonator are particularly 
interesting. Significant research has recently been per- 
formed in this direction (see e.g Ref . [f| @, 0, HI) , by using 
a coupling between the SQUID dynamics and the res- 
onator's mechanical vibrations due to the constraint set 
by the flux quantization phenomenon. Here we will con- 
sider the new possibility for NEM coupling that occurs 
if the nanomechanical element is an integral part of the 
superconducting weak link. In this case the NEM vibra- 
tions directly affect the Cooper pair tunneling and signif- 
icantly modify the properties of the link. With a voltage 



FIG. 1: Sketch of system considered. A nanowire is suspended 
between two superconducting electrodes separated by a trench 
of length L. When the system is biased by a voltage V, the 
Lorentz force caused by the coupling of the Josephson current 
and a transverse magnetic field, H , induces wire vibrations 
described by the coordinate u(x,t). The nonlinear coupling 
leads to a multistability of the system resulting in different dc 
Josephson current regimes (see text). 



The Hamiltonian describing the electronic subsystem 
in the specific model system shown in Fig. [Dreads, 
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h 2 ( d ieHu(x,t)\ 

Ha = A(x) {a x cos </>(£) + sign{x)<7 y sin <fr(t)) , 

where *f?^(x) [^f(x)] are two-component Nambu-spinors 
and <7j are the Pauli matrices in Nambu space [l3]. The 
deflection of the tube is given by u(x,t) = u(x)a(t), 
where u{x) is the normalized (dimensionless) profile of 
the fundamental bending mode and a(t) determines its 
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amplitude; other modes are less important and will be ig- 
nored. The potential U (x) describes the barrier between 
the nanowire and the bulk superconducting electrodes, 
where the gap parameter is A(x) = Ao#(2|x| — L) with 
Ao ~10meV. The phase difference across the junction 
due to the bias voltage V is (j>{t) — 2eVt/h. 

A convenient gauge transformation, see 0] for a similar 
analysis, shifts the vector potential induced by the nan- 
otube deflection from the kinetic part of the Hamiltonian 
to the phase difference between the leads, so that 4>(t) —> 
ip(t) — <f>(t) — a(t)4eH J Q u(x)dx/H. In the adiabatic 
limit, hDip(t) <§; Ao, with D the transparency of the 
barriers, one can then evaluate the fixed-phase ground 
state energy of the electronic subsystem as E(ip) — 
-A [l - Dsm 2 {(p/2)] 1 / 2 and find that the force ex- 
erted on the wire, F — —dE(<p(a))/da, is proportional to 
the Josephson current j = (2e/h)dE(<p)/d<p. The result- 
ing effective equation of motion for the nanowire vibrat- 
ing in its fundamental bending mode describes a forced 
nonlinear oscillator with damping. In terms of the di- 
mensionless coordinate Y(t) = [4eLH/h]a(t) one finds in 
the low transparency limit, D<1, the result 



and 2 (as well as a small peak at V — 1/2), where the 



Y + jY + Y = e sin(v?) 
ip = V -Y . 



(2a) 
(2b) 



Here, 7 = 7/7710; is a dimensionless damping coefficient, 
while e = 8eL 2 H 2 j c / (mHui 2 ) is the amplitude and ip the 
frequency of the driving force with V = 2eV/(Hu>), lo the 
mechanical eigenfrequency, m the mass of the nanowire, 
j c = DAoe/(2h) the critical current and time t mea- 
sured in units of 1/lu. In (|2a[) . the driving force on the 
nanowire is naturally interpreted as the Lorentz force due 
to the coupling between the Josephson current and the 
magnetic field. According to (|2bj) . the phase difference ip 
between the leads evolves in time under the influence of 
both the bias voltage and the electromotive force induced 
by the motion of the wire in the static magnetic field. 

Multiplying (|2aj) with Y and averaging over time we 
find (using the definition of the Josephson current above) 
that in the stationary regime the dc current through the 
system is 
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where (...) denotes time-averaged quantities. To proceed 
with our analysis we consider the specific case of a single- 
wall carbon nanotube wire of diameter 1 nm suspended 
over a length L ~ 1 /im. With j c ~ 100 nA \tj§ one then 
finds that e ~ 3xl0 _3 in a magnetic field of H ~ 20 mT. 
Since 7 = 1/Q, where the quality factor Q ~ 1000, both 
parameters may be considered small, e, 7 -C 1. 

Numerical simulations of the nanowire dynamics using 
the equation of motion ((2]) with Y(0) = Y(0) = show 
distinct resonance peaks in the vibration amplitude at 
integer values of V. Figure [H e.g., shows peaks at V—l 



onset of the V=2 peak depends on the ratio e/7 (14 1. 
For small vibration amplitudes, when one can expand 
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FIG. 2: Time average of the rms nanowire deflection coordi- 
nate Y(t) from a numerical simulation of Eq. Q as a function 
of bias voltage V for different force parameters e (7 = 0.001). 



s'm(Vt + Y) to linear order in Y, these results can read- 
ily be attributed to a direct resonance at V = 1 and a 
parametric resonance atV — 2. In this limit there is a re- 
semblance between the resonances in our system and the 
familiar Fiske effect in Josephson junctions coupled to an 
electromagnetic resonator 12j. However, as can be seen 



in Fig. 2(b) the oscillation amplitude is too large for the 
linear approximation to hold if the driving force is large, 
e > 7. As will be shown below, the resonances in this 
nonlinear regime are significantly different from those of 
the Fiske effect and demonstrate a variety of unusual pe- 
culiarities which could be useful for device applications. 

To analyze the nonlinear regime in the vicinity of the 
resonance peaks it is convenient to apply perturbation 
theory and expand in the small parameters e and 7. With 
this in mind, it is useful to write the dimensionless de- 
flection coordinate Y(t) of the nanowire as 



Y(t) = y/ljt) sin — + 
V n 



Vt , X n(t) 



(4) 



where the amplitude A n {t) = \fl n (t) and phase Xn(t) 
vary slowly in time; I n (t) , Xn(t) ~ 7, e <C 1. Substituting 
this Ansatz into ^ and integrating over the fast oscilla- 
tions one gets two coupled equations for /„ (t) and \ n (t) , 



In = -jl n - 2enJ n (yJTn) sin Xn (5a) 
Xn = S - 2enJ4(V%0 cos Xn ■ (5b) 

Here, J„ are Bessel functions of order n, J' n (\fL^) = 
dJ n {\fln} / dl n and S — V — n. Stationary nonlinear os- 
cillation regimes can now be found by studying the sta- 
tionary points of ([5|) in terms of the system parameters. 

We start our analysis by considering the case of exact 
resonance, (5 = 0, for which Eq. guarantees that a 
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stationary solution given by /„ = always exists. How- 
ever, since J n {x) ~ x n for small x, one immediately finds 
that this solution is unstable for n — 1 (resonance exci- 
tation) , while for n = 2 (parametric excitation) it is only 
unstable if e > 2j. This turns out to be the main differ- 
ence between the two resonance types; the corresponding 
finite-amplitude stationary regimes are qualitatively very 
similar. In the following analysis we will therefore focus 
on the parametric resonance at n = 2, and omit the index 
n on amplitudes, phases and Bessel functions. 

From Eq. {5bJ it is evident that exactly on resonance 
a finite-amplitude stationary regime may be realized ei- 
ther by fixing the phase, cos% = 0, or the amplitude, 
J'(A) = 0. We will refer to these different regimes as 
type I and type II. From (|5a|) it follows that a type-I 
stationary point exists for any e > t\ = 2j. The os- 
cillation amplitude is implicitly given by the equation 
7/e = AJ (A\) / A 2 , which always has a solution in the 
relevant range of parameters, see Fig. 3(a) [13]. Further- 
more, type-II stationary points corresponding to fixed- 
amplitude oscillations, An = A®, where J'(A ) — 0, only 
exists if e > en = ; jAq/AJ(Ao). In this case there exists 
two stationary points of equal amplitude, Afj = Aq, but 
different phases, Xn = 37r/2± arccos(eu/e) 




(a) (b) 

FIG. 3: (a) Plots for solution of Eq. flSaJ with t n = for 
Xn = 3n/2 (type-I regime) and e = 0.005 (red solid), e = 0.003 
(black dashed) and e = 0.001 (green dashed dotted). Cross- 
ings with the straight line correspond to stationary points, 
(b) Numerical solution on resonance of Eq. (|2ap for the time- 
averaged rms nanowire deflection coordinate as a function of 
e. Inset shows corresponding plot for dc current as a function 
of magnetic field. 5 = and 7 = 0.001 throughout. 



A stability analysis shows that the type-I stationary 
point is stable if e < en, but unstable (a saddle point) 
for e > en- The type-II stationary points, on the other 
hand, are always stable if they exist, i.e. when e > en. 
This means that if one increases e, by turning up the 
magnetic field, the nanotubc vibration amplitude will be 
zero (to an accuracy of order 5,7) until e ~ e\. As e is 
varied from e\ to en the amplitude increases from to 
Aq, where it saturates as we increase the magnetic field 
further. This analysis, which also explains the onset of 



vibration amplitude at varying values of e, as shown in 
Fig. 3(b) The inset shows the dc current as a function 
of magnetic field, where H\\ is defined from en oc H 2 ^ as 
above. As the dc current scales as jdc oc (Y(t) 2 }/H 2 one 
finds that the current initially grows with increasing mag- 
netic field strength, pumping energy into the nanoscale 
vibrations, but falls off as 1/H 2 once H > H\\ and the 
vibration amplitude has saturated at Aq. 

Moving off the resonance, 5 becomes non-zero and if 
e > en the degeneracy of the amplitudes A^ at the type- 
II stationary points is lifted. If 5 > the amplitude 
A^j(5) is larger and A^S) smaller than the on-resonance 
value Aq, as shown in Fig. [4] while if S < the op- 
posite is true. As the degeneracy is lifted, the stable 
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the second peak in Fig. 2(b) has been fully confirmed 



by numerically solving the equation of motion @ for the 



FIG. 4: Phase space diagrams in the type-II regime with e = 
0.01 and 7 = 0.001 showing (a) two asymptotically stable 
points with vibration amplitude Ao ~ 3 [dashed line] and (b) 
their shifts from Ao off resonance. 



type-II stationary point that moves to higher amplitudes 
merges with the type-I saddle point and disappears at 
some critical value ±5 C . Consequently, in the interval 
{—S c , 5 C ) there are two different stable nonlinear regimes 
(±) characterized by different nanotube oscillation am- 
plitudes and as a consequence by different dc currents 
through the system. A detailed analysis shows that if 
e — en << en the width 2S C of this window of bistabil- 
ity is cx (e — en) 3 / 2 , while the maximum difference in 
amplitudes \A^ (±S C ) — A^ I (±8 C )\ is cx (e — en). 

The stationary point that describes the system in a 
particular situation depends on the initial conditions. If 
initially A « 0, the system always moves to the stationary 
point with lowest amplitude as the parametric resonance 
develops, i.e. A^ if 5 > and ^ if 5 < 0. However, if 
the system starts from inside the separatrix defining the 
higher than on-resonance stationary point (see Fig. [4]), it 
will achieve a stationary amplitude that is larger than on 
resonance. Alternatively, the system can reach this point 
if the voltage is slowly changed from resonance, as the 
system will follow the trajectory of the stationary point 
at which it is defined exactly on resonance. This repre- 
sents a unique sensitivity in our system to small changes 
in the applied bias voltage. Since the dc current through 
the system depends on the vibration amplitude, jd c oc A 2 , 
it follows that we can predict a non-single valued I — V 
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curve close to resonance. The result is a hysteretic be- 
havior, the origin of which lies in the multistability of the 
pumped nanomechanical vibrations. This means that the 
magnitude of the dc Josephson current in our device is 
sensitive to the pumping history. Such memory effects 
may be employed for different device applications where 
the sensitivity of the nanomechanical initial conditions 
and the possibility to switch the system between two sta- 
ble regimes of vibration can be employed for both sensing 
and memory devices. As an example we discuss briefly 
below how a memory device could work. 
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FIG. 5: Diagram illustrating how small voltage pulses switch 
the dc Josephson current. The two pulses in the inset have 
different effects. After pulse (ii) (blue) the current is the same 
as before, while pulse (i) (red) switches the current by a mea- 
surable amount (from point 1 to 2) thereby storing one bit of 
retrievable information in the device (e = 0.05, 7 = 0.001). 



A scheme for the electrical manipulation of our super- 
conducting nanovibrator is presented in Fig.[5j where the 
starting position 1 corresponds to a bias voltage which is 
slightly off resonance and a nanowire that oscillates with 
an amplitude A^{5 > 0) smaller than on resonance, see 
Fig. |4(b)| Now consider the effect of the voltage pulses 
(i) and (ii) shown in the inset. Pulse (i) moves the sys- 
tem along trajectory (i) to where the vibration ampli- 
tude A n (5 < 0) is larger than on resonance. However, 
at 5 = — <5 C this asymptotically stable point merges with 
the third stationary point (saddle) and becomes unstable. 
The system therefore jumps to the second asymptotically 
stable point, where the vibration amplitude Af^S < 0) is 
smaller than on resonance. When the voltage is increased 
again, the system will move to position 2, where the vi- 
bration amplitude Aj^(S > 0) and hence the dc current 
is larger than at position 1. One concludes that pulse (i) 
writes one bit of information, which is stored as a mea- 
surably larger dc current. Pulse (ii), on the other hand, 
moves the initial stability point back and forth along tra- 
jectory (ii) and returns it to the initial position 1. 

It is interesting to again compare the phenomena dis- 



cussed in this Letter with the Fiske effect 12[. Repeat- 
ing our analysis we find the low-amplitude behavior to 
be similar for the two systems. However, we predict that 
a dynamical multistability will appear at a certain value, 
en, of the driving Lorentz force. This does not occur in 
the Fiske effect, where the vibration amplitude exactly 
on resonance follows the stable solution corresponding to 
cos x = in (|5b|) for all driving forces e. 

To conclude we have shown that for a nanowire sus- 
pended between two voltage-biased superconducting elec- 
trodes in a transverse magnetic field, pronounced reso- 
nance phenomena can be found at discrete values of the 
driving voltage. Our analysis shows that the behavior of 
this system is governed by an effective equation of mo- 
tion whose solution gives the amplitude of the nanowire 
oscillations and the dc Josephson current as a function 
of system parameters. Most importantly, it was shown 
that for realistic experimental parameters the system can 
be driven into a multistable regime by varying the mag- 
netic field strength. The possibility to pump energy into 
the mechanical vibrations of a suspended nanowire and 
the ensuing dynamical multistability of the vibration am- 
plitude and dc current makes this superconducting nano- 
electromechanical device a unique system, where the sen- 
sitivity to initial conditions and switching between two 
stable regimes can be probed experimentally. 
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